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Abstract: The Gauss law needs to be imposed on quantum states to guarantee gauge 
invariance when one studies gauge theory in hamiltonian formalism. In this work, we 
propose an efficient variational method based on the matrix product ansatz for a Zi lattice 
gauge theory on a spatial ladder chain. Gauge invariant low-lying states are identified by 
evaluating expectation values of the Gauss law operator after numerical diagonalization of 
the gauge hamiltonian. 
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1. Introduction 

The importance of the first-principle study in quantum chromodynamics is increasing 
largely because RHIC experiment has started and LHC is also coming. For precise de- 
scription of high-energy heavy ion collisions, gauge theory needs to be studied at finite 
temperature and density in a systematic way. Ideally, we should also have a methodology 
for tracing time-evolution of quantum states based on the Schrddinger equation because 
heavy ion collisions should be treated as non-equilibrium evolving systems rather than 
static. Lattice gauge theory is the most useful method for studying the quark-gluon sys- 
tems at zero and finite temperature. However, Monte Carlo integration does not work 
for lattice gauge theory with large chemical potential because of the severe sign prob- 
lem. It would be worthwhile to pursue a systematic variational approach to gauge theory 
|3], [|, J5|. In this paper, we apply the matrix product ansatz to a Z2 hamiltonian lattice 
gauge theory on a spatial ladder lattice based on the previous work in a U(l) lattice gauge 
theory ||. 

The matrix product ansatz [[?J is a variational method that originates in DMRG (den- 



sity matrix renormalization group) |8[ ||, 10] . DMRG has been developed as the method that 
gives the most accurate results for spin and fermion chain models such as one-dimensional 
quantum Heisenberg and Hubbard models at zero and finite temperature jll|, 12, 13]. 1 



DMRG is also useful for diagonalization of transfer matrices in two-dimensional classical 
statistical systems [14|. DMRG has been extended to two-dimensional quantum systems 



15, [R]] and can work for bosonic degrees of freedom p7| ]. Since DMRG is a variational 
method based on diagonalization of hamiltonian and transfer matrices, it is free from the 
sign problem. Actually, DMRG has been successful in accurate study of ^-vacuum in the 



By "d-dimensional" , we mean (1 + d)-dimensional spacetime. 
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massive Schwinger model [18|. There is an old prediction by S. Coleman that quarks decon- 
fine at 9 = tt [19]. However, the model has not been analyzed accurately with the Monte 
Carlo method because the topological terms give complex action in the Euclidean theory. 

The matrix product ansatz is a result of large simplification of trial wavefunction based 
on the knowledge established in the past DMRG studies. Therefore, the ansatz takes over 
the good points of DMRG. Besides, the matrix product ansatz has advantage to DMRG 
because the former can treat periodic one-dimensional systems accurately [^]. Recent 
interesting progress of the matrix product ansatz is its application to non-equilibrium 
quantum physics and quantum information theory [21]. The matrix product ansatz is a 
promising approach to further refinement and extension of the past DMRG studies. 

In general, calculation of expectation values of hamiltonian becomes difficult exponen- 
tially as the system size increases. If the matrix product ansatz is introduced, the energy 
function has a simple matrix product form, which can be evaluated easily using a com- 
puter if the matrix size is small. It is expected that exponentially difficult problems are 
reduced into small tractable ones by the matrix product ansatz. Actually, the ansatz has 
been successful in giving accurate results in many one-dimensional quantum systems, where 
calculation errors can be controlled systematically. We can say that the matrix product 
ansatz is a first-principle variational method. 

Lattice gauge hamiltonian is obtained by choosing temporal gauge in partition function 
of Euclidean lattice gauge theory [22|. In hamiltonian formalism, gauge invariance needs 
to be maintained explicitly by imposing the Gauss law on the Hilbert space. It is a 
hard task to construct gauge-singlet variational space for general gauge group [p3[| . On 
the other hand, Euclidean lattice gauge theory can keep gauge invariance manifestly by 
construction. This is one of the reasons why hamiltonian version of lattice gauge theory 
is not popular. In addition, no systematic methods had been known for diagonalization of 
gauge hamiltonian before the matrix product ansatz was applied to lattice gauge theory 
in ref. ||. If trial wavefunction is constrained directly with the Gauss law, the advantage 
of the matrix product ansatz is completely spoiled because calculation of energy function 
becomes impossible in a practical sense. If the hamiltonian is diagonalized without the 
Gauss law, all possible states are obtained including gauge variant states. However, it 
must be possible to extract gauge invariant states because all eigenstates of the hamiltonian 
can be classified using generators of the considered gauge group. Therefore, if the matrix 
product ansatz is used, we better start from the whole Hilbert space and then identify 
gauge invariant states using the Gauss law operator after all calculations. 

In hamiltonian lattice gauge theory, the dimension of spatial lattice needs to be two or 
larger because the plaquette operator is a two-dimensional object. In this paper, we study 
a Z2 lattice gauge theory on a spatial ladder chain for simplicity. The length of the chain 
needs to be sufficiently long for the matrix product ansatz to work well. The ladder chain 
is squashed into a chain because one-dimensional structure needs to be found in order to 
use the matrix product ansatz. 

In the previous work || , we have studied one- and two-dimensional 5 = 1/2 Heisenberg 
models and a U(l) lattice gauge theory on a ladder chain using the matrix product ansatz 
in a similar way, where energy function is minimized using the Powell method. When the 
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number of parameters is very large, such naive minimization is not useful because it takes 
long time to reach the bottom of the energy function. In this work, we use a diagonalization 



method introduced in ref. [20 to obtain sufficient accuracy. 

This paper is organized as follows. In section ^, hamiltonian lattice formulation of 
the Z<i lattice gauge theory is briefly reviewed. In section [3|, the matrix product ansatz is 
introduced and applied to a Z2 lattice gauge hamiltonian. In the original construction, the 
matrix product states is assumed to have translational invariance. In this work, that con- 
dition is not imposed on variational space before diagonalization of the gauge hamiltonian. 
In section ||], numerical results are given. Section [5] is devoted to summary. 

2. Quantum hamiltonian in the Z 2 lattice gauge theory 

We are going to introduce the Z2 lattice gauge theory, which was invented by F. Wegner 



As seen in the literature, the simplicity of the model is useful for testing a new idea [25, 26]. 
The model cannot have non- vanishing magnetization because local gauge symmetry cannot 
break spontaneously, which is known as the Elitzur's theorem [^7]]. However, the model 
can have nontrivial phases depending on dimensionality. 

We are interested in quantum hamiltonian of the model. Statistical mechanics and 
quantum hamiltonian are connected through the transfer matrix formalism. In the Z2 
lattice gauge theory, quantum hamiltonian is obtained by choosing temporal gauge in the 
partition function [^| 

# = -5>x(M-A^P(n,i,j), (2.1) 

n,i n,i,j 

where a x and a z are spin operators 

= (0 1\ = (1 

°* V 1 0/ ° z V° - 1 

and P is a plaquette operator 

P(n,i,j) = a z (n,i)a z (n + i,j)a z (n + i + j,-i)a z (n + (2.2) 

In eq. (|2.1|), the first and second summations are taken on the spatial lattice for all possible 
link and plaquette operators, respectively. In general, arbitrary states can be represented 
as a superposition of products of |±) n ,i> which are eigenstates of the spin operator a z (n,i) 

<jg(n,i)\±) n>i = ±|±) n> i. 

Let us introduce time-independent operators G(n), each of which flips spins on all the 
links emerging from a site n 

G(n) = l[a x (n,i). (2.3) 

±i 

We have 

G^ 1 {n)a x (m, i)G(n) = a x (m,i), 
G _1 (n)cr z (n,i)G(n) = -a z (n,i), 
G^ 1 (n)a z (m,i)G(n) = a z (m,i), 
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where the last formula applies only if the link (m, i) is not contained in G(n). The operator 
G(n) defines local gauge transformation 

G{n)- l HG{n) = H. (2.4) 

In order for physical quantities to be gauge invariant, quantum states need to be invariant 
under gauge transformation 

G(n)|tf) = I*). (2.5) 

We need to impose the Gauss law Q2.5| ) on the wavefunction to keep gauge invariance. 
Otherwise, unphysical states may be obtained because gauge invariance is not guaranteed. 
When a state \*S>) satisfies the Gauss law (|2.5|), magnetization vanishes because a relation 
(9\a g (n)\-9) = -(¥|<7*(n)|tt) holds. 



3. Matrix product ansatz on a ladder lattice 

We are going to introduce the matrix product ansatz, which is a variational method in- 
spired from density matrix renormalization group (DMRG). DMRG is a variational method 



that can reproduces very accurate results in one-dimensional quantum systems ||, ||, 17]. 
In DMRG, wavefunction is represented as a product of orthogonal matrices because ba- 
sis states are rotated for optimization with orthogonal matrices that diagonalize density 
matrices. The success of DMRG allows us to parametrize wavefunction as a product of 
finite-dimensional matrices from the beginning. This simplification of wavefunction is called 
the matrix product ansatz Although DMRG has slow convergence in one-dimensional 
quantum systems with periodic boundary conditions, the matrix product ansatz gives much 
better accuracy [20 ]. 



A[ Sl ] A[s 2 ] A[s L ] A[ Sl ] 

Figure 1: A spatial chain with lattice size L. The open circles indicate periodicity. The site 
variables are dynamical. The same set of matrices A[s) is assigned to all the sites. 

In the original construction, the matrix product state is parametrized as follows Q: 

\n=l s„ / 

where {|sn)} is a complete set of basis states for the n-th site (see figure [J). Hamiltonian 
needs to have periodicity for consistency with the trace operation in the variational state. 
As a result, energy becomes a function of the matrices A[s]. The minimum of the energy 
function corresponds to the ground state. We are going to apply the ansatz to a Z<i gauge 
theory and see its compatibility with gauge symmetry. 

Since this work is the first application of the matrix product ansatz to Z<i gauge theory, 
we would like to consider a simple model. The simplest one is a Z2 hamiltonian lattice gauge 
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theory on a spatial ladder lattice (see figure |2|). We assume periodicity in the horizontal 
direction on the ladder for later convenience. In figure ^, periodicity is denoted with the 
open circles. 




C L [u L ] 



Figure 2: A spatial ladder chain with lattice size L. The open circles indicate periodicity. The 
link variables are dynamical. Different sets of matrices are assigned to links. 



In the free case A = 0, the hamiltonian fl2.1|) can be diagonalized analytically. The 
vacuum state is given by 

JJ-L(|+) nt< + |-) nii ). (3.2) 



vac) 



Vacuum expectation values of the hamiltonian and plaquette operators are (H) /L = —3 
and (P) = 0, respectively. In this case, one-dimensional matrices are sufficient to represent 
the vacuum state (|3.2| ) 

A[+] = A{-] = -L (3-3) 



When the coupling constant is very large A — > oo, the first term of the hamiltonian 
( |2.1| ) can be neglected. As a result, the hamiltonian has a diagonal form. The vacuum 
state is given by 



vac) 



1 2L 

m=0 n.i 



(3.4) 



where the operator G^ represents a summation of all possible products that are composed 
of m pieces of different Gaussdaw operators G. The operator G^ 71 ' has 2LC m terms (G^ = 
1). In eq. ( |3.4[ ), the two definitions are identical. Then, we have (H)/L — > —A and (P) — > 1 
in the limit A — > oo. The states (^2) and (3^) are both gauge invariant, G(n)|vac) = |vac). 

The Z2 lattice gauge model has only link variables. In our construction, each link is 
assigned a different set of matrices A n , B n , and C n for parametrization of wavefunction 
(see figure §) . The index n represents the n-th square on the ladder chain and runs from 
1 to L. The dimension of the matrices is M. Then, our matrix product state is 



|*)=tr [J!] Yl E A n[sn]B n [t n ]C n [u n ] \s n ) n \t 

\n=l s n =± t n =± u n =± 



n/n\ Un |n 1 



(3.5) 



where the matrices are multiplied in ascending order keeping the order of A n B n C n , and 
the basis states |s) n , \t) n , and \u) n are eigenstates of the spin operator u z as before. In this 
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expression, the variables s, t, and u are used instead of the index % to denote the position 
of the links. The implementation of the matrix product ansatz means that a ladder lattice 
has been represented as a one-dimensional system with non-nearest neighbor interactions. 
Gauge invariance of matrix product states will be discussed in the next section. 

If we require orthogonality of optimum basis states according to ref . |7]] , we have 

M 

£ X)(Xn[s])tf(*n[*])i'j = *«', (3-6) 
3=1 s=± 
M 

Y,Y,( X ^M x n[s]k> = 5jf, (3-7) 

i=l s=± 

where X stands for A, B, and C. If these conditions are not imposed, norm of the matrix 
product state ( |3.5| ) may becomes very small, which results in numerical instability. 
Energy 

E= V ' . ' 7 , 3.8 
(*|^) 

is a function of the matrices ^4 n [s], 5 n [t], and C n [u]. The numerator and denominator can 
be calculated by evaluating trace of a product of 3L matrices numerically: 

L / n— 1+L mod L 

($?\H\V) = - tr ( a nynZ n + X n b n Z n + X n y n C n ) Y[ W m + 

n=l \ m=n+l mod L 

n— 1+L mod L \ 



+ \ot. n f3 ri r y n Ot n +\ mod LUn+1 mod L z n+1 mod L ] | Wm ] , (3.9) 

m=n+2 mod L / 

=tr m u; m ^ , (3.10) 

where 

"n = ^foW^M ® A n [s'], j3 n = ^2{a z ) tt 'B* n {t] ® B n [t'], 7n = ^(^WC*M C n [it / ], 

s,s' t,t' u,u' 



t 

w n = x n y n z n . 



The dimension of the matrices on the left hand side is M 2 . By the outer product symbol 
(g>, we mean 

a (i,k),(j,i) = y~](7ss'Aij[s\A k i[s']. 



The minimum of the energy function ( |3.8| ) corresponds to the ground state, which can be 
obtained based on matrix diagonalization as explained below. We can reduce the mini- 
mization problem ( |3.8[ ) into a generalized eigenvalue problem [20| 

v^Hv = Ev^Nv, (3.11) 
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where H and N are 2M 2 by 2M 2 matrices. To understand what is going here, let us 
consider how energy can be minimized by varying A n [s] when other matrices are fixed. 
Note that equations ( |3^9| ) and ( |3.10|) are bilinear of the matrix ^4 n [s] 

<*|#|*> = Y, T / ( A ni S ])iAid,sUkU)( A n[t\)kl, (3-12) 
i,j,k,l s,t 

<*i*> = y E«W)«^(vi.-).(*A*)(^w)«' ( 3 - 13 ) 

i,j,k,l s 

where the matrix N is diagonal for the indices s and t. Once these expressions are obtained 
and the variational parameters A n [s] are regarded as a vector v, the minimization problem 
|OD reduces to ( gjjj ). 

There is one more trick that needs to be implemented. As explained in equations 
Q3.6| ) and Q3.7| ), we encounter numerical instability on the right hand side in eq. ( |3.11| ) if 
the above procedure is used as it is. This is because the matrix N may have very small 
eigenvalues if the matrices A n [s], -B n [s], and C n [s] are varied freely. For this reason, we 
need to impose one of the conditions ( |3.6| ) or ( |3.7| ) on the matrices. 

A matrix can be decomposed into a product of three matrices, which is called singu- 
lar value decomposition. For example, if we regard a tensor A[s] as a matrix A^u )g \ = 
(A[s])ij = Au, singular value decomposition of Ajj is given by 

M 

Au = Y U IK D K V KJ . (3.14) 

K=l 

where Dk are the singular values of the matrix A pS|]. The matrices U and V are or- 
thonormal 

M 2M 

Uj K U IK > = 8 KK i, Y V KJ V K'J = &KK'- (3.15) 

1=1 J=l 

The decomposition A = U'V, where U' = UD is a square matrix, is the key of the trick. 
Consider a part of the matrix product wavefunction 

C[u]A[s] = C'[u]A'[s], (3.16) 



where (A'[s])jj = Vi u s \ and C'[u] = C[u]U'. If the decomposition ( 0.14 ) is accurate, the 



both representations give the same result. Based on this trick, the eigenvalue problem 



(3.11) is solved successively for all the sets of the matrices starting from the right end in 
figure § 

Cl[ul] - B L [t L ] - A L [s L ] > dim] -> Btih] Ai[si], (3.17) 



which we call sweep. In one sweep process, eq. ( 3.11f) is solved 3L times. If this sweep 



process is repeated several times, energy of low-lying states converges to some value. Cal- 
culation error can be controlled systematically by increasing the matrix dimension M. See 
Appendix |A| for generation of initial matrices. 
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4. Numerical results 



The matrix product ansatz assumes large lattice. Our lattice size L = 500 is sufficiently 
large. We solve the generalized eigenvalue problem ( p. 11 ) using LAPACK |3(J. For steady 



states, real matrices are sufficient for parameterizing the matrix product state ( |3.5| ). Con- 
vergence of energy needs to be checked for the number of sweeps and the matrix dimension 
M. Energy density E/L converges in accuracy of five digits or higher after two sweeps 
when the matrix size M is fixed. 



M 


E /L 


Ex/L 


E 2 /L 


E 3 /L 


E A /L 


E 5 /L 


A = 0.1 


2 


-3.001 


-2.997 


-2.997 


-2.997 


-2.993 


-2.993 


3 


-3.001 


-2.997 


-2.997 


-2.997 


-2.994 


-2.993 


4 


-3.001 


-2.997 


-2.997 


-2.997 


-2.997 


-2.995 


A = 1 


2 


-3.124 


-3.121 


-3.121 


-3.118 


-3.114 


-3.112 


3 


-3.124 


-3.121 


-3.121 


-3.118 


-3.114 


-3.112 


4 


-3.124 


-3.121 


-3.121 


-3.118 


-3.114 


-3.112 


A = 10 


2 


-10.27 


-10.27 


-10.27 


-10.27 


-10.23 


-10.23 


3 


-10.27 


-10.27 


-10.27 


-10.27 


-10.26 


-10.23 


4 


-10.27 


-10.27 


-10.27 


-10.27 


-10.26 


-10.23 



Table 1: Energy density E/L of six low-lying states is listed for A = 0.1, 1, and 10 when lattice 
size is L = 500. Good convergence of energy is obtained with small M. 



Table [T] shows energy spectra of six low-lying states for three values of the coupling 
constant: A = 0.1, 1, and 10. The sweep process has been repeated twice. In this model, 
convergence of energy is very fast in contrast to Heisenberg chains HQ- Small matrix 
dimension is sufficient for good convergence. Since we have obtained low-lying states with- 
out imposing the Gauss law on the variational space, gauge variant states are contained. 
In table |l[ gauge invariant states are denoted with underlines. The other states are gauge 
variant. As we will see, gauge invariant physical states can be identified by calculating 
expectation values of the Gauss law operator. 

In the ladder chain model, the Gauss law operator G{n) is a product of three a z 
operators (two horizontal and one vertical). We evaluate expectation values of G(n) on the 
upper lattice sites shown in figure ^. Then, the number of the Gauss law operators to be 
evaluated is L. Expectation values on the lower sites are same as the upper ones because of 
reflection symmetry. Figures [3| plots expectation values of the Gauss law operator (G(n)) 
in the case of A = 10 for the states (a) Eq, (b) E\, (c) E2, and (d) E3. In figures |3] (a) and 
(d), the Gauss law G(n) = 1 is satisfied uniformly on every lattice sites. Therefore, the 
obtained states Eq and E3 are gauge invariant. On the other hand, in figures || (b) and (c), 
the states E\ and E2 are gauge variant because gauge symmetry is definitely broken at the 
site n = 500. The position of this special lattice site depends on where the sweep process 
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(a) The state E 

1 

(G) o 
-1 

100 200 300 400 500 

n 

(b) The state E\ 

1 

(G) o 
-1 

100 200 300 400 500 

n 

(c) The state E 2 

i — 



100 200 300 400 500 

n 

(d) The state E 3 



1 

(G) o 
-1 

100 200 300 400 500 

n 

Figure 3: Expectation values of the Gauss law operator is plotted for the four low-ling states (a) 
E Q , (b) Ei, (c) E 2 , and (d) E 3 with X = 10, L = 500, and M = 4. The circles are the calculated 
values. The states (a) and (d) are gauge invariant because the Gauss law is satisfied on every lattice 
sites. On the other hand, the states (b) and (c) are gauge variant because (G(500)) = — 1. These 
statements hold in accuracy of seven digits or higher. 

ends. The relation (G(n)) = 1 or —1 holds for the obtained low-lying states in accuracy of 
seven digits or higher when M = 4. 

According to the Elizur's theorem, gauge variant operators have vanishing expectation 
values. We have checked that expectation values of single spin operators cr z (n) vanish for 
the gauge invariant states in accuracy of ten digits or higher when M = 4. We also have 
checked that the gauge variant states have vanishing expectation values of cr z (n) in the 
same accuracy. These statements apply to the low-lying states shown in table [j]. 

In this way, we can classify the obtained states into gauge invariant states and others. 
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E /L 




2 4 6 8 10 



A 

Figure 4: Vacuum energy density Eq/L is plotted as a function of the coupling constant A with 
lattice size L = 500. The circles are the calculated values, which are consistent with the exact ones 
given in sec. ||. 

Figure || plots vacuum energy density as a function of the coupling constant A. 

In order to see how spatial distribution look like, we evaluate expectation values of the 
plaquette operator for the vacuum and the excited states. Figure || plots expectation values 
of the plaquette as functions of the spatial lattice coordinate n for the gauge invariant low- 
lying states (a) Eq, (b) £3, and (c) E$ in the case of A = 10. As expected, the vacuum 
state Eq has complete uniformity. On the other hand, the excited states E3 and E§ have 
lumps around n = 500. The higher excited states have similar lumps around the boundary. 
The obtained solutions have periodicity because the plaquette distributions are continuous 
around the boundary. The position of the lumps can be moved without changing energy 
in the sweep process as explained in figure ||. Figure ||| plots vacuum expectation values of 
the plaquette operator as a function of the coupling constant A. For small and large A, the 
tendency of energy and plaquette is consistent with the exact values given in section ||. 

5. Summary 

We have extracted gauge invariant physical states in a Z2 hamiltonian lattice gauge theory 
on a spatial ladder chain. The calculations are based on the matrix product ansatz, which 
gives sufficiently convergent energy and wavefunction for low- lying states. In the future 
studies, similar calculations should be tested in higher dimensional lattice gauge theory 
including supersymmetric cases The proposed method will be useful especially for 

non-perturbative analysis of vacuum structure at the amplitude level. 
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(a) The vacuum state E 



(P) 



(P) 



0.5 





(P) 0.5 








100 200 300 

n 


400 


500 


(b) The excited state E 3 








o 



I I I I 


i 




o 


100 200 300 

n 


400 


500 


(c) The excited state E 5 






i i i i i i 


o 




o 


o 




o 



100 



200 300 

n 



400 



500 



Figure 5: Vacuum expectation values of the plaquette operator are plotted for the gauge invariant 
states (a) Eq, (b) £3, and (c) £5 with A = 10 and L = 500. The excited states (b) and (c) have 
non-uniformity around the boundary n = 500. 

A. Generation of initial matrices 

Before starting sweep process, we need to prepare initial values of matrices for the vari- 
ational state (|3.5|). As explained in section p| the matrices need to satisfy one of the 



orthogonality conditions ( |3.6| ) or (3.7) for numerical stability. The algorithm shown below 
is also useful for minimization of energy function with the Powell method ||. 

Let us consider real matrices A[s] that satisfy the following normalization condition 



M K 

j = l 8=1 



(A.l) 



where M is the dimension of the matrices A[s] and K is the degrees of freedom of each site 
or link. To parametrize the matrices A[s], we introduce .fTM-dimensional vectors 



(n) 



(n) 

' • • • ' V KM 



_„ +1 ,0,...,0), n 



1,...,M. 



(A.2) 



n—l 



These vectors are linearly independent and can be orthonormalized using the Gram- 
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A 

Figure 6: Vacuum expectation value of the plaquette operator is plotted as a function of the 
coupling constant A with lattice size L=500. The circles are the calculated values, which are 
consistent with the exact ones given in sec. |[ 



Schmidt method 



,.(1) 



l«Wl 



6 (*+i) = V (k+D _ ^( v (fc+i) ja W) (») 



n=l 



,(*+!) 



\ b (k+l)\ 



(A.3) 
(A.4) 
(A.5) 



where k = 1, . . . , M — 1 and the brackets represent inner product. The orthonormalized 
vectors are used to parametrize the matrices A[s]jj = Aju s -\ 



(i> 



A; 



*,(?» 



fa 



(A.6) 



which satisfy the conditions (A.l). When eq. ( |A.1| ) is satisfied, the number of independent 
degrees of freedom associated with the matrices A[s] is given by 



KM Z 



M(M + 1) 



which is equal to that of the vectors v ^ 

M-l 

(KM -n)-M. 



(A.l) 



(Ai 



n=0 



In eq. ( A. 8 ), the first term counts the number of the parameters vf^ and the second term 
comes from the normalization conditions ( |A.3| ) and 
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